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I.  INTRODUCTION 


The  conjugate  gradient  method  was  developed  by  Hestenes  and  inde¬ 
pendently  by  Steifel,  but  the  initial  comprehensive  treatment  was  a  joint 
publication(ref.  1).  The  first  application  of  the  method  to  the  synthesis  of 
antenna  patterns  appears  to  have  been  implemented  by  Fong  and  Birgen- 
heier(ref.  2)  to  compute  the  line  source  required  to  produce  a  desired  power 
pattern.  This  application  concerned  the  solution  of  a  Fredholm  equation  of 
the  first  kind  by  minimizing  an  appropriate  weighted  least-squares  integral. 
The  algorithm  used  in  reference  2  is  a  computationally  fast  nonlinear  min¬ 
imization  scheme,  but  it  has  two  major  disadvantages.  First,  it  may  not 
converge  for  many  functions  that  represent  the  power  pattern  of  a  discrete 
array.  Second,  an  external  one-dimensional  minimization  algorithm  must  be 
supplied  by  the  user.  These  two  disadvantages  can  lead  to  an  unpredictable 
breakdown  in  the  algorithm. 

This  study  has  two  major  parts.  The  first  part  is  the  development  of  a 
reliable  way  to  apply  a  conjugate  gradient  algorithm  to  the  minimization  of 
an  arbitrary  nonlinear  function.  The  second  part  is  the  application  of  that 
algorithm  to  the  minimization  of  the  maximum  sidelobe  level  of  a  planar 
array  of  dipoles  using  phase-only  synthesis. 
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II.  GEOMETRY  AND  FAR  FIELD 

The  array  configuration  chosen  for  this  study  is  a  two-dimensional  array 
of  half-wavelength  horizontal  dipoles  placed  above  a  ground  plane.  The 
dipoles  are  placed  such  that  there  is  quadrant  symmetry,  and  a  circular 
perimeter  is  approximated.  The  top  view  of  a  sample  array  is  shown  in  Fig. 
1. 


Fig.  1.  Geometry  of  Planar  Array 


Assuming  an  eJwf  time  convention,  the  total  far-zone  electric  field  due  to 
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a  single  ^-directed  current  source  is  given  by 

-jk0R 


E  =  i  *°^iT|sin 


in(0)^  4-  cos(0)  cos(<jf>)^J  Nx 


(2.1) 


where 


and 


Nx  =  ///  -M*.  y>  z)ej(',*+vv+'/k,°-u,"’z)dxdydz  (2.2) 


u  =  ko  sin(fi)  cos(0) 
v  =  ko  sin(0)  sin(0). 


(2.3) 

(2.4) 


If  the  source  consists  of  4 N  equal-amplitude,  half- wavelength  dipoles  placed 
a  distance  h  above  a  ground  plane  and  if  the  kih  current  has  the  form 


=  h6{y  -  yh)S(z) sin  k0  -  \x  -  xk 

then  equation  (2.2)  reduces  to 

Nx  =  j4k0  sin  (hy/kl-u*-v?)  g  hej 


I*  -  aefc |  <  y  (2.5) 


«»(£:») 

The  far-zone  power  pattern  is  then  given  by 


l7k+Vn(«i”)] 


(2.6) 


(\ N  \2  /4 N  \: 

Y  h  cos[7*  +  i>h{u,  v)]  )  +  sin [7*  +  rpk{u,  w)]  ) 
U=i  /  Vk=i  / 


P[y,u,v)  =  f3(u,v) 

where 

0{U'v)  =  klh^S'm2  (^o-«2-^)  cos2 


(2.7) 


(2.8) 


and  i’k(u,v)  =  uxk  +  rj/^.  Because  quadrant  symmetry  is  assumed,  the  power 
pattern  simplifies  to 


P(7,u,u)  =  0(u,v) 
where  (k  -  zos(uxk)  cos(vyk). 


\*  =  1 


)  COS  +  ( y^/'‘C*(u-i’)si'>(7ik) 


^=1 


(2.9) 


G 


Antenna  pattern  synthesis  involves  controlling  the  power  pattern.  For  a 
known  amplitude  excitation,  control  of  the  power  pattern  may  be  achieved 
by  adjusting  the  relative  phase  excitation  of  each  element.  The  goal  of  the 
synthesis  will  be  to  maximize  the  power  in  the  main  beam  while  simulta¬ 
neously  minimizing  the  peak  sidelobe  power.  Observation  of  equation  (2.9) 
indicates  that  the  far-zone  power  pattern  is  a  nonlinear  function  of  the  un¬ 
known  excitation  phase.  Therefore,  a  nonlinear  minimization  scheme  will  be 
required  in  the  synthesis. 


III.  A  NONLINEAR  CONJUGATE  GRADIENT  METHOD 


As  the  previous  section  indicated,  minimizing  the  sidelobe  level  power 
of  an  antenna  pattern  involves  the  minimization  of  a  nonlinear  function.  A 
conjugate  gradient  method  may  be  used  to  find  the  minimum  value  of  the 
nonlinear  function  /O),  wrhere 

7  =  [7i  72  1n}T  (3.1) 

with  each  7*  representing  the  relative  phase  shift  for  the  fcth  element.  It  will 
be  assumed  that  /( 7)  is  a  nonconvex,  twice-difierentiable  real  function  of  a 
real  variable.  This  implies  that  there  are  multiple  solutions  to  the  equation 

V/  =  0  (3.2) 

where  V/  is  the  gradient  of  /( 7)  defined  by 

*/=[&  ■  ^\T ■  (3-3) 

Each  solution  of  equation  (3.2)  is  called  a  “stationary”  or  “critical”  point. 
A  conjugate  gradient  solution  is  possible  by  assuming  an  initial  estimate  71 
and  approximating  /( 7)  as  a  quadratic  truncation  of  a  Taylor  series  expanded 
around  this  point  as 

/( 7)  =  /(71)  -  (7  -  7 l)TG  +  ^(7  -  7')rtf(7  -  71)  +  0(17  -  T1 13)  (3.4) 

where  G  is  the  negative  gradient  of  /( 7)  evaluated  at  71  given  by 

c-  -(W)(r')=  If^a1)  &(,•)  &b'>r  (5.5) 


and  //  i 

s  the 

Hessian  of  /( 

7)  evaluated 

at  71  given  by 

r 

-Al/_  (7i|l 

II 

(vr/ih1 1 

n'J  ,,■» 

1  1  ' 

nil,  (7‘  )  ■  • 

-?V  (TI ) 

(375 

. 1,3 1  (V  ) 

>±‘l  '  | 

(I 


It  is  important  to  note  that  the  Hessian  matrix  H  is  always  symmetric  but  not 
necessarily  positive  definite.  The  matrix  would  be  positive  definite  if,  given 
any  nonzero  vector  p,  the  condition  { Hp,p )  >  0  is  satisfied.  The  brackets 
denote  an  inner  product  between  two  N  x  1  vectors  defined  by 

N 

(3  7) 

n  =  1 

Making  the  substitution  z  =  7  -71  and  excluding  the  constant  and  higher 
order  terms  of  equation  (3.4)  yields  the  function 

c(z)  =  -ztHz  ztC,  (3.8) 

which,  when  minimized,  produces  the  same  solution  as  solving  the  linear 
system  of  equations  Hz  —  G.  Because  H  and  G  are  implicit  functions  of  z,  the 
algorithm  must  be  restarted  periodically  to  yield  the  correct  solution  to  the 
original  problem.  If  //  is  positive  definite,  then  a  possible  algorithm  which 
iterates  K  times  with  L  restarts  is  given  as  follows: 


begin  loop  for  l=l,...,L 


compute  Gi  =  -(V/)(yl),  Hi  =  und  set  21  ~  0 


begin  loop  for  k  K 


r 1  -  G, 

(3.9) 

do  ~  1  .  , \ 

(r'.r1) 

(3.10) 

P1  Ho?' 

(3.11) 

<** 

(HiPh  P*) 

(.*•  1  - 

4  1 

If  .  k 

1 3.13 

-  Z  1  Uj <,p 

/41 

Tk  CtkHlpk 

Cut 

/  0  1  r 

err 

(r* , rk) 

(0.1  .» 
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if  err  <  tolerance,  termii  ate  k  loop,  else 


1 


(3.16) 

(3.17) 


Pk 

k  t  1 


(r*  +  i,rk  +  l) 
Pk  t  Pkrk  +  l 


continue  k  loop 


7 


i  +  i 


7l  4 


,*4i 


(3.18) 


continue  l  loop 

This  algorithm  requires  that  //  he  explicitly  known  in  order  to  compute  a* 
and  rku.  An  alternate  algorithm  that  does  not  require  computation  of  H 
may  be  derived  by  noting  that  rk  +  l  -  -(Vc)(:fc)  «  -(V  f)(ylc ).  Also,  a k  may  be 
approximated  by  solving  the  one-dimensional  minimization  problem 

choose  a*  to  minimize  /( 7*  f  a*p*)  (3.19) 

This  yields  an  alternate  algorithm  proposed  in  references  2  and  3  which  do~s 
not  require  the  explicit  calculation  of  H .  The  algorithm  is  given  as  follows: 


begin  loop  for  l  1,..  ,L 


compute  Ci  —  (V/)( 71),  Hi  —  (VV/)( 7f)  and  set  z1  =  0 

r1  =  G, 

(rbr1) 
p'  --  do r' 


begin  loop  for  k  / . A 


rl.nns r  minimize  /(->*  1  ) 


h 


i 


i 


k 


i  ^  /)(->'  i  .-A") 

.*11  r*  1 1 . 

,  r* i 


(3.20) 

(3.21) 

(3.22) 


(3.231 


(3.21) 


(3. 2S! 


(3.21.) 


M 


if  err  <  tol'  .ncc,  terminate  k  loop,  else 


1 


(3.27) 

(3.28) 


'  (r*  +  >,r*  +  >) 

pk  +  l  =  pk  +  0::rk  +  1 

continue  k  loop 

7'+l  =  7'  4  zk  +  1  (3.29) 

continue  l  loop 

This  algorithm  avoids  computation  of  the  Hessian  but  may  not  always 
yield  a  solution.  If  H  is  not  positive  definite,  then  the  scalar  ak  may  be  infi¬ 
nite  or  negative,  which  violates  the  conditions  under  which  the  algorithm  is 
derived.  This  may  cause  an  unpredictable  breakdown  in  the  one-dimensional 
minimization  algorithm  used  to  find  equation  (3.23).  Thus,  although  the  algo¬ 
rithm  presented  in  equations  (3.20) -  (3.29)  is  computationally  efficient,  the  low 
reliability  would  require  a  constant  user  interface  and  an  additional  uecision 
concerning  the  contingency  of  breakdown  of  the  algorithm.  An  alternative 
algorithm  that  does  not  suffer  from  these  diadvantages  is  presented  next. 

If  the  Hessian  is  nonsingular,  but  possibly  indefinite,  then  an  algorithm 
that  iterates  K  times  with  L  restarts  is  given  as  follows: 

begin  loop  for  l-l,...,L 


compute  Gi  -  —(V /)( 71 ),  Hi  ~  (VV f){ 7*)  and  set  z1  —  0 

r1  -  Ct'i 

/in  -  ,  -  — 1  -  - 

,  H 1  r 1  ,  H 1  r 1 ) 

p1  linHir] 


(3.30) 

(3.31) 

(3.32) 


In  gin  loop  for  h  I .  .  h 


'n  {Hirk.  Hip'  ) 


(3.33) 
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2*+1  =zk  +  akpk 

(3.34) 

r^‘  -r*-Qfc7/,p* 

(3.35) 

(r‘  +  »,r‘  +  l) 
err  ~  - - — — — - 

(rk,rk) 

(3.36) 

if  err  <  tolerance,  terminate  k  loop,  else 

(3.37) 

p*  +  l  +  1 

(3.38) 

continue  k  loop 

7,+  1  ^7!  +  2*  +  ! 

(3.39) 

continue  l  loop 

This  algorithm  requires  knowledge  of  both  the  gradient  and  Hessian  of  the 
function  but  offers  greater  reliability  against  failure. 

The  quadratic  is  a  local  approximation,  so  that  if  the  initial  guess  is  not 
near  a  solution  to  the  global  problem,  then  there  is  a  point  of  diminishing  re¬ 
turns  in  terms  of  minimizing  the  quadratic.  In  other  words,  if  the  quadratic 
is  not  near  a  minimum  of  the  original  function,  then  it  is  not  beneficial 
to  minimize  it  very  far.  It  is  better  to  perform  a  small  number  of  conju¬ 
gate  gradient  iterations  on  the  quadratic  and  then  restart.  As  the  quadratic 
approximation  gets  closer  to  the  actual  trough  of  the  original  function,  it  be¬ 
comes  advantageous  to  increase  the  number  of  conjugate  gradient  iterations. 
In  general,  the  number  of  restarts  L  will  be  a  large  number  such  that  L  >>  !. 
This  contributes  to  the  majority  of  the  computation,  because  it  takes  much 
more  time  to  compute  the  gradient  and  Hessian  than  to  perform  a  conjugate 
gradient  iteration  Tin-  conservative  approach  would  be  to  always  keep  l\ 
mall  such  t hat  l  •  A  lit.  However,  the  optimal  approach  would  bo  to  allow 
A  to  be  a  function  of  /,  ,  such  that  A  A  ( /.)  could  bo  incremented  when  the 
quadratic  gets  closer  to  a  solution  region. 
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IV.  MINIMIZING  MAXIMUM  SIDELOBE  LEVEL 


Now  that  the  numerical  algorithm  has  been  established,  it  is  necessary  to 
find  the  appropriate  function  to  minimize.  The  desired  goal  of  the  synthesis 
is  to  minimize  the  peak  sidelobe  power  level  while  keeping  the  peak  of  the 
main  beam  as  high  as  possible.  It  is  well  known  that  the  minimum  peak 
sidelobe  level  for  a  given  beam  width  occurs  when  all  the  sidelobes  are  at  the 
same  height.  Therefore,  in  order  to  implement  the  algorithm  of  the  previous 
section,  it  is  necessary  to  find  a  suitable  function  that,  when  minimized,  will 
yield  this  behavior.  The  function  used  by  Deford  and  Gandhi  (refs.  4  and 
5)  is  given  by 


/( 7)  = 


P(j,  v.,  r.) 

>(7,0,0) 


where  P(-y,Ui,Vi)  is  the  peak  sidelobe  power  level  located  at  the  point  ( uitVi ) 
and  P(y,  0, 0)  is  the  power  at  the  main  beam  peak.  If  this  function  is  used, 
then  as  the  iterations  proceed,  a  new  peak  sidelobe  is  found  and  the  function 
is  changed.  Unfortunately,  using  this  approach  can  lead  to  stability  problems 
as  the  sidelobe  level  reaches  a  uniform  level.  This  is  because  as  the  sidelobe 


level  becomes  more  uniform,  the  pattern  becomes  very  sensitive  to  a  change 
in  any  sidelobe.  This  can  lead  to  oscillations  and  instability. 

The  approach  proposed  in  this  report  is  to  minimize  the  weighted  average 
of  the  current  peak  sidelobe  level  with  all  previous  peak  sidelobe  levels  at 
each  restart  of  the  algorithm  given  by  equations  (3.30)  through  (3.39).  More 
specifically,  at  each  restart,  the  sidelobe  region  is  searched  to  find  the  peak 
sidelobe  point.  If  the  peak  sidelobe  point  occurs  at  a  previous  peak,  then 
the  components  of  the  weighted  average  remain  unchanged.  However,  if  the 
p<’nk  occurs  at  a  m*w  sample  point,  then  a  new  function  is  formed  bv  taking 
lie-  weighted  average  of  the  new  point  with  all  the  previous  peak  points. 
As  an  example,  suppose  that  at  some  restart  number,  the  sidelobe  region  is 
searched  and  the  point  (it*,?,*),  k  /  i  is  found  to  have  the  highest  sidelobe. 


If) 


Then  equation 


Iteration  1 


Iteration 


Iteration 


Fig.  2.  Weighted  Average  Minimization 
4.1)  would  be  modified  to 


/( 7)  = 


ELi 


P(7,0,0)  J22l=l 


p2{ 

P{ 7,  «i,  Vi) 


(4.2) 
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The  averaging  effect  is  chosen  because  as  the  sidelobes  approach  a  uniform 
level,  the  averaging  provides  an  automatic  damping  effect  that  reduces  the 
sensitivity  of  the  sidelobe  pattern  to  small  changes  in  the  phase.  The  aver¬ 
age  is  weighted  so  that  more  emphasis  is  placed  on  minimizing  the  largest 
sidelobes. 

Figure  2  exhibits  the  typical  evolution  of  an  axial  cut  in  the  sidelobe 
pattern,  assuming  that  the  synthesis  produced  uniform  sidelobes  in  five  iter¬ 
ations.  At  the  first  iteration,  the  near-in  sidelobe  has  the  highest  peak.  After 
the  first  minimization,  the  same  point  still  has  the  highest  peak.  At  the  third 
iteration,  the  second  sidelobe  now  has  the  highest  level,  so  a  weighted  aver¬ 
age  of  the  two  lobes  is  formed.  At  the  fourth  iteration,  the  fourth  lobe  now 
is  the  highest,  and  finally  at  the  fifth  iteration,  all  the  lobes  are  at  an  equal 
level.  Because  of  the  conservation  of  energy,  this  level  will  always  be  between 
the  highest  and  lowest  peaks  of  the  pattern  with  uniform  illumination. 

A  general  weighted  average  function  may  now  be  defined  as 


/( 7) 


=  Pg(  T) 

Po(l) 


where  P„( 7)  is  the  weighted  average  sidelobe  level  at  the  Ith  restart 
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and  the  u  -  v  plane  is  sampled  as  in  Fig.  3.  The  components  of  the 
and  Hessian  are  computed  from 
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Fig.  3.  Sampled  u  -  v  Plane 


the  first  and  second  derivatives  may  be  expanded  as 
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and  letting  ippq  =  uqxp  +  vqyp  yields 
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Note  that  for  the  broadside  array  with  a  main  beam  at  9  =  0  and  <j>  —  0 
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dim  dln 

where  CP?  =  cos)^*,)  cos^y,). 

The  phase  may  be  initialized  as  the  sum  of  a  linear  shift  and  a.  pertur¬ 
bation  term  given  by 


Ik  -  2jt (/*,)**  t  /yo yfc )  t  £(**,y*)  (4.22) 

when-  siiilfl,,) r»s((/>„)  and  /yf,  siii(fln) sin(</>0).  The  term  £(**,?/*)  must 

be  added  if  all  amplitudes  are  equal  because  observation  of  equation  (4.20) 
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indicates  that  a  linear  phase  is  a  stationary  point  of  equation  (4.3).  Deford 
and  Gandhi  (ref.  4)  have  found  that  although  the  final  phase  distribution  is 
strongly  dependent  on  the  initial  guess,  the  resulting  peak  sidelobe  level  is 
not . 

It  is  important  from  a  practical  standpoint  to  have  a  measure  of  the 
ability  of  the  array  to  concentrate  energy  in  the  main  beam.  The  parameter 
of  interest  is  the  directivity  defined  by 


Directivity  = 


^^ATEPTE 


(4.23) 


where  Aa  is  the  aperture  area,  ATE  is  the  amplitude  taper  efficiency,  and  PTE 
is  the  phase  taper  efficiency.  The  ATE  is  a  measure  of  how  the  amplitude 
variations  of  the  electric  field  across  the  aperture  of  the  array  affect  the 
concentration  of  energy  into  the  main  beam.  Letting  »  denote  the  surface, 
the  ATE  is  given  by 


Aaff.\E\>d* 


(eLia)2 

"ZLxH ' 


(4.24) 


The  PTE  is  a  measure  of  how  the  phase  variations  of  the  electric  field  across 
the  aperture  of  the  array  affect  the  concentration  of  energy  into  the  main 
beam.  The  PTE  is  given  by 


pTE_.  i  If.Ed'l2 

iff,  \E\2ds)2 

[ELi  I*  cos(t* )]  +  [£"=1  Ik  sin(7*  )] 


(4.25) 


A  phase-only  array  with  a  uniform  amplitude  distribution  has  an  ATE  1. 
This  implies  that,  the  maximum  power  in  the  main  beam  is  achieved  by  a 
uniform  amplitude  and  phase  distribution. 
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A  common  characteristic  of  both  amplitude  and  phase  synthesis  is  that 
as  the  sidelobe  level  is  lowered,  the  beam  width  increases.  This  increase  in 
beam  width  may  be  measured  with  respect  to  the  beam  width  of  a  uniform 
array.  Assuming  a  circular  aperture  of  diameter  D,  if  D  >>  A0,  then  the  half 
power  beam  width  (HPBW)  of  a  uniform  array  is  given  by 

HPBW  =  58.95°^  (4.26) 

and  the  beam  width  between  first  nulls  (BWFN)  is  given  by 

BWFN  =  139.76°  —  .  (4.27) 

Note  that  the  beam  width  is  primarily  determined  by  the  diameter  of  the 
array.  It  can  be  observed  by  Fourier  transform  theory  that  the  HPBW 
and  the  BWFN  for  a  uniform  amplitude  and  phase  array  are  the  minimum 
achievable. 
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V.  RESULTS 


All  the  results  shown  in  this  report  were  computed  by  the  Fortran  77 
program  CGPHSA.  A  description  of  the  input  parameters  as  well  as  the  source 
code  listing  is  given  in  the  appendix.  All  results  were  computed  on  an  IBM 
3090  computer. 

Figure  4  shows  the  peak  sidelobe  level  as  a  function  of  the  array  diameter 
for  a  uniform  phase  array  and  a  phase- synthesized  array.  The  dipole  length 
is  given  by  <fj  =  0.5Ao,  and  the  center- to- center  spacing  is  given  by  a  =  0.577Ao. 
The  dipoles  are  placed  at  a  distance  of  h  =  0.25Ao  above  a  perfectly  conducting 
ground  plane.  The  number  of  restarts  was  approximately  L  —  ION,  where  N 
is  the  number  of  elements  in  the  first  quadrant.  The  number  of  conjugate 
gradient  iterations  was  set  at  K  —  2.  The  sidelobe  level  of  the  uniform 
phase  array  approaches  the  limit  of  the  continuous  distribution,  which  is 
approximately  -17.6  dB.  The  synthesized  phase  array  sidelobe  level  indicates 
a  logarithmic  dependence  on  the  array  diameter.  A  logarithmic  curve  fit 
yields  the  equation 

SLL(dB)  =  -4.831n(D)  -  10.22.  (5.1) 

Extrapolation  based  on  this  equation  indicates  that  a  sidelobe  level  of  -40 
dB  would  require  about  a  470Ao  diameter  aperture. 

The  relative  aperture  efficiency  is  shown  in  Fig.  5.  The  curve  follows 
the  same  logarithmic  behavior  as  the  peak  sidelobe.  At  D  =  20Ao,  the  power 
in  the  main  beam  has  dropped  to  about  48%  of  what  it  is  for  a  uniform 
distribution.  This  level  is  similar  to  that  of  an  equivalent  array  using  only 
amplitude  tapering. 

Figures  6  and  7  compare  the  half  power  beamwidth  (HPBW7)  and 
beaimvidth  between  first  nulls  (BWFN)  of  the  synthesized  array  to  a  uniform 
array.  The  synthesized  curves  appear  to  be  just-scaled  versions  of  equations 
(  1.26)  and  (4.27).  The  HPBW  for  the  synthesized  array  appears  to  be  about. 
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Diameter  (wavelengths) 


Fig.  6.  Half-Power  Beam  Width 
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Fig.  7.  Beam  Width  Between  First  Nulls 
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1.08  times  the  uniform  HPBW,  and  the  BWFN  for  the  synthesized  array 
is  about  1.1  times  that  of  the  uniform  phase  case.  This  beam  broadening 
is  also  a  characteristic  of  amplitude  tapered  arrays  and  appears  to  be  at  a 
comparable  level. 

Figures  8  and  9  are  power  patterns  of  a  D  —  7A0  diameter  array  plotted 
in  decibels  in  cylindrical  coordinates  on  a  circular  base  with  a  -32  dB  floor. 
Figures  10  and  11  are  the  same  patterns  plotted  in  spherical  coordinates. 
As  shown  in  Figs.  8  and  10,  the  sidelobe  level  of  a  uniform  phase  array 
decreases  monotonically  as  0  goes  from  0°  to  90°.  Figures  9  and  11  indicate 
that  the  lowest  peak  sidelobe  level  occurs  when  all  the  sidelobes  have  the 
same  height.  Therefore,  the  synthesis  of  low  sidelobes  results  in  a  uniform 
distribution  of  energy  in  all  sidelobes.  The  total  energy  is  conserved,  because 
the  energy  in  the  peak  sidelobes  and  the  main  beam  is  distributed  over  the 
lowest  sidelobes  to  bri  >g  all  the  sidelobe  peaks  to  an  equal  height. 


Fig.  8.  Uniform  Phase  Power  Pattern  in  Cylindrical  Coordinates 
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Fig.  10.  Uniform  Phase  Power  Pattern  in  Spherical  Coordinates 
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Fig.  11.  Synthesized  Phase  Power  Pattern  in  Spherical  Coordinates 
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VI.  CONCLUSIONS 


The  application  of  the  conjugate  gradient  method  based  algorithm  de¬ 
scribed  in  this  report  appears  to  yield  a  very  reliable  numerical  procedure. 
If  an  accurate  initial  guess  is  not  available,  then  the  number  of  restarts  must 
be  large  and  the  number  of  conjugate  gradient  iterations  per  restart  must 
be  small.  The  results  indicate  that  the  peak  sidelobe  power  level  may  be 
decreased  by  phase-only  synthesis.  The  level  appears  to  be  a  logarithmic 
function  of  the  array  diameter,  which  limits  the  usefulness  to  very  large  ar¬ 
rays.  The  convergence  is  stabilized  by  forming  a  weighted  average  of  present 
and  past  peak  sidelobe  power  points. 
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APPENDIX:  FORTRAN  77  PROGRAM 


INPUTS 

The  Fortran  77  program  CGPHSA  synthesizes  the  phase  required  to  yield  the 
lowest  sidelobe  level  for  a  circular  two-dimensional  array  of  x-directed  half¬ 
wavelength  dipoles.  Because  Fortran  77  allows  only  static  storage  allocation, 
the  dimensions  of  all  vectors  must  appear  in  the  calling  program.  The  required 

dimension  parameters  and  declarations  are  given  as  follows. 

1  c**********************************************************************c 

2  PROGRAM  CGPHSA 

3  PARAMETER  (M=146,NS=980) 

4  REAL*4  X(M),T(M),A(M),P(M),W(M),R(M),Q(M),Z(M)IG(M),E(M),T(M) 

5  A.UCNSj.VdfSj.BtNS) ,H(M,M) ,S(M,NS) ,C(M,M) 

6  INTEGER*4  J(NS) 

7  D=7 . 0 

8  SI= . 677 

9  SY= .  6 

10  CALL  GEOMET(D,SX,SY,N,X,T) 

11  N0R=0 

12  CALL  PARBOL(D,NOR,NfX,Y,A) 

13  CALL  CGPATKM.N.X.T.D.A.ItS^.W.R.Q.Z.G.U.V.T.B.J.E.H.S.C) 

14 

15 
18 

17  END 

18  C********4t*****+***************+*****************+**+****+****+*+******C 
The  maximum  vector  dimension  M  is  the  number  of  elements  that  fits  into 

a  square  aperture  of  side  D ,  multiplied  by  the  ratio  of  the  areas  of  a  circle  of 
diameter  D  to  a  square  of  side  D,  and  is  given  by 

M  =  nearest  integer  — - -  (.4.1) 

The  parameter  NT  should  be  set  to 

NT  -  nearest  integer  [6. 4411  -  5]  (.4.2) 

and  tin-  parameter  NS  is  then  set  to 

NS *NT(9  I  NT)  (,4.;$) 
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The  subroutine  GEOMET  returns  the  number  of  elements  and  rectangular  coor¬ 
dinates  of  the  center  point  of  all  dipoles  in  the  first  quadrant  of  the  array.  The 
diameter  P,  center-to-center  x-spacing  SI,  and  the  center-to-center  y-spacing  ST 
are  the  required  inputs.  The  following  inputs  are  passed  to  the  subroutine  GEOMET. 


Input 

Type 

Description 

D 

R4 

diameter  of  the  aperture  in  wavelengths 

SX 

R4 

center-to-center  dipole  spacing  in  the  x  direction 

ST 

R4 

center-to-center  dipole  spacing  in  the  y  direction 

N 

14 

unspecified  integer 

X(M) 

R4 

unspecified  vector 

T(M) 

R4 

unspecified  vector 

The  following  output: 

s  are  returned  from  the  subroutine  GEOMET. 

Output 

Type 

Description 

N 

14 

number  of  dipoles  in  the  first  quadrant 

X(M) 

R4 

x  coordinates  of  the  center  of  each  dipole 

T(M) 

R4 

y  coordinates  of  the  center  of  each  dipole 

The  subroutine  parbol  returns  the  parabolic  amplitude  distribution  of  the 

elements. 

The  diameter  and  the  order  of  the  parabolic  taper  NOR  are  the  required 

inputs.  The  following  inputs  are  passed  to  the  subroutine  PARBOL 

Input 

Type 

Description 

D 

R4 

diameter  of  the  aperture  in  wavelengths 

NOR 

14 

order  of  the  parabolic  taper 

N 

14 

number  of  dipoles  in  the  first  quadrant 

KM) 

R4 

x  coordinates  of  the  center  of  each  dipole 

T(M) 

R4 

y  coordinates  of  the  center  of  each  dipole 

A(M) 

R4 

unspecified  vector 

The  following  output 

is  returned  from  the  subroutine  PARBOL. 

Output 

Typ- 

Description 

A  ( M ) 

R4 

current  amplitudes  of  each  dipole 
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The  following  inputs  are  passed  to  the  subroutine  cgpatt. 


Input 

Type 

Description 

D 

R4 

diameter  of  the  aperture  in  wavelengths 

M 

14 

maximum  dimension  of  vectors 

N 

14 

number  of  dipoles  in  the  first  quadrant 

NS 

14 

total  number  of  sample  points  in  the  u-v  plane 

X(M) 

R4 

x  coordinates  of  the  center  of  each  dipole 

Y(M) 

R4 

y  coordinates  of  the  center  of  each  dipole 

A(M) 

R4 

current  amplitudes  of  each  dipole 

P(M) 

R4 

scratch  vector 

W(M) 

R4 

scratch  vector 

R(M) 

R4 

scratch  vector 

Q(M) 

R4 

scratch  vector 

Z(M) 

R4 

scratch  vector 

G(M) 

R4 

scratch  vector 

T(M) 

R4 

scratch  vector 

U(NS) 

R4 

scratch  vector 

V(NS) 

R4 

scratch  vector 

J(NS) 

R4 

scratch  vector 

B(NS) 

R4 

scratch  vector 

E(M) 

R4 

scratch  matrix 

R4 

scratch  matrix 

S(M,NS) 

R4 

scratch  matrix 

C(M.M) 

R4 

^cratch  matrix 

OUTPUTS 


The  folowing  outputs  are  returned  from  the  subroutine  CGPATT. 


Output 

Description 

P(N) 

relative  phases  of  each  dipole 

T(N) 

relative  phases  of  each  dipole 

Z(l) 

half-power  beam  width 

Z(2) 

beam  width  between  first  nulls 

Z(3) 

amplitude  taper  efficiency 

Z(4) 

phase  taper  efficiency 

Unit 


radians 

degrees 

degrees 

degrees 

percent 

percent 
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LISTING 


1  PROGRAM  CGPHSA 

2  C*********************************+'‘*******++*****************+*+******C 

3  C  THIS  PROGRAM  COMPUTES  THE  PHASES  REQUIRED  FOR  A  CIRCULAR  PLAHAR  C 

4  C  ARRAY  OF  X-DIRECTED  HALF  WAVELENGTH  DIPOLES  PLACED  A  QUARTER  OF  A  C 

5  C  WAVELENGTH  ABOVE  A  GROUND  PLANE.  C 

e  c*******************************#*************************************^ 

7  C  TIMOTHY  J.  PETERS  LAST  UPDATED  C 

s  C  THE  AEROSPACE  CORPORATION  OS/12/89  C 

g  C  ELECTRONICS  RESEARCH  LABORATORY  C 

to  C  ELECTROMAGNETIC  SCIENCES  DEPARTMENT  C 

u  C  2350  EAST  EL  SEGUNDO  BOULEVARD  C 

12  C  EL  SEGUNDO,  CA  90245  C 

13  C*******  +  iM*iMiM**4^4^******4^*4^4^*  +  *  +  **v''-t**************************+C 

14  C  VARIABLE  DESCRIPTIONS:  C 


IB 

c 

C 

10 

c 

NAME 

TYPE 

DESCRIPTION 

C 

17 

c 

c 

18 

c 

D 

R4 

ARRAY  DIAMETER  IN  WAVELENGTHS 

c 

19 

c 

SX 

R4 

CENTER  TO  CENTER  DIPOLE  SPACING  IN  THE  X  DIRECTION 

c 

20 

c 

SY 

R4 

CENTER  TO  CENTER  DIPOLE  SPACING  IN  THE  Y  DIRECTION 

c 

21 

c 

M 

14 

MAXIMUM  DIMENSION  OF  VECTORS  COMPUTED  FROM  THE  FORMULA 

c 

22 

c 

M=NEAREST  INTEGER(PI*D*D/(16*SX*SY) ) 

c 

23 

c 

N 

14 

NUMBER  OF  UNKNOWN  ELEMENTS  (N  <=  M) 

c 

24 

c 

NT 

14 

NUMBER  OF  RADIAL  SAMPLE  RINGS  COMPUTED  FRDM  THE  FORMULA 

c 

26 

c 

NT=NEAREST  INTEGER ( 8 . 44*D-5 ) 

c 

20 

c 

MO 

14 

NUMBER  OF  SAMPLE  POINTS  IN  U-V  PLANE  COMPUTED  FROM 

c 

27 

c 

FORMULA  NS=NT(9+NT)/2 

c 

28 

c 

X(M) 

R4 

I  COORDINATES  OF  THE  CENTER  OF  EACH  DIPOLE 

c 

29 

c 

Y(M) 

R4 

Y  COORDINATE  OF  THE  CENTER  OF  EACH  DIPOLE 

c 

30 

C  A(M) 

R4 

AMPLITUDES  OF  EACH  DIPOLE 

c 

31 

C 

P(M) 

R4 

SCRATCH  VECTOR 

c 

32 

c 

W(M) 

R4 

SCRATCH  VECTOR 

c 

33 

c 

R(M) 

R4 

SCRATCH  VECTOR 

c 

34 

c  q(M) 

R4 

SCRATCH  VECTOR 

c 

35 

c 

Z(M) 

R4 

SCRATCH  VECTOR 

c 

30 

c 

G(M) 

R4 

SCRATCH  VECTOR 

c 

37 

c 

T(M) 

R4 

SCRATCH  VECTOR 

c 

38 

c  U(NS) 

R4 

SCRATCH  VECTOR 

c 

39 

c 

V(NS ) 

R4 

SCRATCH  VECTOR 

c 

40 

c 

J(NS) 

14 

SCRATCH  VECTOR 

c 

1 1 

c 

B  ( NS ) 

R4 

SCRATCH  VECTOR 

c 

i  :* 

c 

E(K) 

R4 

SCRATCH  VECTOR 

c 

4  3 

n 

H(M,M) 

R4 

SCRATCH  MATRIX 

c 

44 

c 

S ( M , NS) 

R4 

SCRATCH  MATRIX 

c 

45 

c 

C(K.M) 

R4 

SCRATCH  MATRIX 

c 

40 


47  PARAMETER  (M=29 ,NS=980) 

48  REAL*4  X(M),Y(M),A(M),P(M) ,W(M) ,R(M) ,Q(M) ,Z(M) ,G(M) ,T(M) ,D(NS) 

49  A,V(NS),B(NS),E(M),H(M,M),S(M,NS),C(M,M) 

so  INTEGER*4  J(NS) 

si  0PEN(UNIT=2 ,FILE= ’CGPOUT’ ,  STATUS = ’UNKNOWN’ ) 

52  C**************************************** ***************************** *C 
63  C  COMPUTE  THE  COORDINATES  OF  EACH  DIPOLE  EXCLUDING  THE  REFERENCE.  C 

54  C  THE  TOTAL  NUMBER  OF  DIPOLES  SHOULD  ALREADY  BE  KNOWN  AS  NSE.  C 

55  C**********************************************************************C 

5«  D=7 . 0 

57  SX=0 . S77 

58  SY=0 . S77 

59  CALL  GEOMET(D , SX , SY , N , I , Y) 

80  c********************************************* ************************ *C 
61  C  COMPUTE  THE  EXCITATION  AMPLITUDE  OF  EACH  DIPOLE  C 

82  C* ****************** ***************************************************C 

63  N0R=0 

64  CALL  PARBOL(D,NOR,N,X,Y,A) 

85  C* *************************** ******* ***********************************C 
66  C  COMPUTE  THE  SOLUTION  C 

87  C**********************************************************************C 
68  CALL  CGPATT(M,N,X,T,D,A,NS,P,W,R,Q,Z,G,U,V,T,B,J,E,H,S,C) 

89  C**** **************************************** **************************0 

70  C  PRINT  THE  PHASE  IN  DEGREES.  C 

71  C* *************************************************************** ******C 


72 

WRITE(2 , 103) 

73 

WRITE(2 , 104) 

(I,P(I) ,1=1, N) 

74 

WRITE(2 , 10B) 

Z(l) 

75 

WRITE(2 , 106) 

Z(2) 

70 

WRITE(2 , 107 ) 

Z(3) 

77 

WRITE(2,108) 

Z(4) 

78  C**********************************************************************c 

79  C  FORMATS  C 

80  C**********************************************************************C 

81  103  F0RMAT(1H1,27X, ’COMPUTED  PHASE  OF  EACH  ELEMENT  IN  RADIANS’,// 

82  A , 61 , ’ ELEMENT ’ , 0X , ’ PHASE ' ,71,' ELEMENT ’ 

83  A, flX, ’PHASE’ ,7X, ’ELEMENT' ,6X, ’PHASE’ ,7X, ’ELEMENT’ 

84  A, flX, ’PHASE’ ,/,6X, ’NUMBER’ ,17X 

85  A, ’NUMBER’ ,17X, ’NUMBER’ ,17X 

86  A, ’NUMBER’ ,/) 

87  104  FORMAT( 6X , 14 , 6X , F7 . 2 , 2X , 14 , 6X , F7 . 2 , 2X , 14 , 6X , F7 . 2 , 2X , 14 , 6X , F7 . 2 ) 

ro  105  F0RMAT(10X, ’HALF  POWER  BEAMWIDTH  (DEGREES)  ’  2X.F6.2) 

80  100  F0RMAT(10X, 'BEAMWIDTH  BETWEEN  FIRST  NULLS  (DEGREES)’  2X.F6.2) 

on  107  F0RMAT(10X, ’AMPLITUDE  TAPER  EFFICIENCY  (PERCENT)  ’  2X.F6.2) 

oi  108  FORMATdOX, ’PHASE  TAPER  EFFICIENCY  (PERCENT)  ’  2X.F6.2) 

o:’  END 

03  C 

01  SUBROUTINE  GEOMET(D , SX , ST ,N , X , Y) 
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95  C*************  +  ****************  +  *****  +  **************************','**  +  +  **C 
99  C  THIS  SUBROUTINE  COMPUTES  THE  COORDINATES  OF  THE  CENTER  OF  EACH  C 

97  C  I  DIRECTED  DIPOLE  FILLING  THE  FIRST  QUADRANT  OF  A  CIRCULAR  APERTURE  C 

98  C  OF  DIAMETER  D  WITH  A  DIPOLE  LENGTH  OF  DL.  THE  CENTER  TO  CENTER  C 

99  C  X  SPACING  IS  SX  AND  THE  CENTER  TO  CENTER  T  SPACING  IS  ST.  C 

100  C**********************************************************************C 

101  REALM  I(*),Y(*) 

102  REALM  DX  ,DT  ,XMIN ,  YMIN  ,IS  ,TS  ,D2  ,DL2  ,  AR,RC 

103  DL=0.6 

104  IMIN=SX/2.0 

105  YMIN=SY/2 .0 

109  DX=SX 

107  DY=ST 

108  DL2=DL/2.0 

109  D2=D/2 . 0 

no  SM=AMIN1(SX,SY) 

111  NS=INT(D/(2.0*SM))M 

112  N=0 

113  DO  1  1=0, NS 

114  YS=YMIN+I*DY 

115  DO  2  J=0 , NS 

119  XS=IMIN+J*DI 

117  AR=(XS+DL2)*(XS+DL2)+TS*YS 

118  RC=DSQRT(AR) 

119  IF(RC  .LE.  D2)  THEN 

120  N=N+1 

121  I(N)=IS 

122  Y(N)=TS 

123  ELSE 

124  END  IF 

125  2  CONTINUE 

129  1  CONTINUE 

127  RETURN 

128  END 

129  C 

130  SUBROUTINE  PARBOL(D , NOR.N ,X , Y , A) 

131  c****************************** ***************************** ***********C 

132  C  THIS  SUBROUTINE  COMPUTES  THE  DIPOLE  AMPLITUDE  EXCITATION  USING  A  C 

133  C  PARABOLIC  DISTRIBUTION.  C 

134  C* I****************************************************************** ***C 

135  REALM  X(  +  )  ,Y(  +  ),A(  +  ) 

130  DO  1  1  =  1, N 

137  R=SQRT( X(I)*X(I)+Y(I)*Y(I)) 

I3H  A(  I  )  -  ( 1 .0-(2.0*R/D)**2)**N0R 

i3t»  I  CONTINUE 

i  \«  RETURN 

in  END 

1  12  C 
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M3  SUBROUTINE  CGPATT(MDIM ,NUE,X ,Y ,D , AI ,NMP ,P , W ,R, Q , Z ,G ,U , V, TV , IP ,B 

144  A,PW,H,PSI,JC) 

145  C**********************************************************************G 


140  C  THIS  SUBROUTINE  IMPLEMENTS  A  CONJUGATE  GRADIENT  METHOD.  C 

147  C*** ****************************************************** *************C 

148  C  VARIABLE  DESCRIPTIONS:  C 

14B  C  C 

iso  C  NAME  TYPE  DESCRIPTION  C 

151  C  C 

152  C  NUE  14  NUMBER  0*  UNKNOWN  ELEMENTS  C 

153  C  I(NUE)  R4  I  COORDINATES  OF  THE  CENTER  OF  EACH  DIPOLE  C 

154  C  Y(NUE)  R4  Y  COORDINATE  OF  THE  CENTER  OF  EACH  DIPOLE  C 

165  C  D  R4  DIAMETER  C 

156  C  HT  R4  HEIGHT  OF  THE  DIPOLE  ARRAY  ABOVE  THE  GROUND  PLANE  Z 

157  C  AIUUE)  R4  AMPLITUDES  OF  EACH  DIPOLE  C 

158  C  NMP  14  NUMBER  OF  MATCH  POINTS  IN  U-V  SPACE  C 

iso  C  NJP  14  MAXIMUM  NUMBER  OF  PEAK  SIDELOBE  LEVEL  POINTS  C 

ieo  C  NT  14  NUMBER  OF  RADIAL  SAMPLE  RINGS  C 

101  C  IRS  14  NUMBER  OF  RESTARTS  C 

102  C  P(NUE)  R4  SEARCH  DIRECTION  C 

163  C  W(NUE)  R4  UNKNOWN  PHASE  AT  EACH  DIPOLE  C 

104  C  R(NUE)  R4  RESIDUAL  C 

105  C  q(NUE)  R4  WORK  VECTOR  FOR  CG  ALGORITHM  C 

106  C  Z(NUE)  R4  ACTUAL  UNKNOWN  VECTOR  FOR  THE  CG  ALGORITHM  C 

107  C  G(NUE)  R4  GRADIENT  C 

168  C  U(NMP)  R4  U  COORDINATES  OF  THE  SAMPLE  SPACE  C 

169  C  V(NMP)  R4  V  COORDINATES  OF  THE  SAMPLE  SPACE  C 

170  C  TV (NUE)  R4  TEMPORARY  VECTOR  THAT  STORES  THE  PHASE  AT  LOWEST  SLL  C 

171  C  IP (NMP)  14  TAG  ARRAY  CONTAINING  ALL  PREVIOUS  PEAK  SLL  LOCATIONS  C 

172  C  B(NMP)  R4  BETA  C 

173  C  PW(NJP)  R4  POWER  AT  EACH  PEAK  SLL  LOCATION  C 

174  C  H(NUE ,NUE)  R4  HESSIAN  C 

175  C  PSI(NUE.NMP)  R4  MULTIPLIER  DUE  TO  ELEMENT  PATTERN  C 

178  C  JC(NJP.NUE)  R4  JACOBIAN  C 

177  C************************************** ******************************* *C 

178  C  NMP  EQUALS  THE  SUM  FROM  K=0  TO  K=NT-1  OF  (IPM+K*IPD)  WHICH  EQUALS  C 

179  C  NMP=NT*(2*IPM+(NT-l)*IPD)/2 


180  C****** ************************ ****************************************0 

181  REAL*4  H(MDIM.MDIM) ,PSI(MDIM , NMP ) , X(* ) , Y(* ) , P( * ) , W( * ) ,R( * ) , Q( * ) 

182  »,Z(*),G(*),B(*),U(*) ,V(*),PW(*),JC(MDIM,MDIK),AI(*) ,TV(*) 

183  INTEG£R*4  IP(*) 

1 04  HT-0 . 25 

185  NJP=MDIM 

i  so  NT-NINT(6.44+D-5) 

is,  irs=10*NUE 

100  SLLMIN=1000 . 0 

iso  RAD= . 174S3293E-01 

ion  PI= .3141S93E401 
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181  TP= . 028318E3E+O1 

192  C**********************************************************************C 

193  C  COMPOTE  THE  AMPLITUDE  TAPER  LOSS.  C 

194  C**********************************************************************C 

ias  CALL  AMPTAP ( HOE ,  AI ,  A  JL  ) 

198  C* ******************************************************************** 

197  C  INITIALIZE  THE  ELEMENT  PHASES  USING  A  LINEaR  PHASE  AND  PERTURBATION.  C 

198  C**********************************************************************C 

199  CALL  INITAL(NUE,D,X,T,W) 

200  c************************* ************************************* ********c 

201  C  FIND  THE  BEAMWIDTH  OF  THE  INITIAL  PATTERN.  C 

202  C* *************************************************************** ******C 

203  CALL  BEAMWDCNUE.I.T.AI.W.HT.HPBW.BWFN) 

204  DO  66  IJ=1,NUE 

205  Z(IJ)=0.0 

208  66  CONTINUE 

207  CALL  BEAMWD(NUE ,X , ? ,  AI , Z ,HT ,HPBWZ ,BWFNZ) 

208  C******************** ************************************************* *C 

209  C  COMPUTE  THE  SAMPLE  POINTS  IN  U-P  SPACE.  C 

210  C**********************************************************************C 

211  TMIN=BWFN/2 . 0 

212  TMAI=88 . 0 

213  PMIN=0 . 0 

214  PMAX=90 . 0 

2'  IPM=6 

216  IPD=1 

217  CALL  ANGLES (TMIN ,TMAX,PMIN ,PMAX ,Ni , IPM , IPD ,U, V ) 

218  CALL  BETASO(HT.BO) 

219  CALL  BETASI(NMP,U,V,HT,B) 

220  CALL  PSIMNI(MDIM ,NUE ,NMP ,X,T,U,TfPSI) 

221  CALL  POWERO (NUE ,Z , AI ,B0 , P01 ) 

222  INIT=0 

423  CALL  POUMAX(MDIM , NUE ,NMP ,  INIT , Z ,PSI ,B , AI , IP ,NP , LMAXi , PMAI1 ) 

224  P1DB=10.0*AL0G10(PMAII/P01) 

225  INIT=0 

228  CALL  POWMAX (MDIE , NUE , NMP , INIT , W , PS I , B , AI , IP , NP , LMAXI , PMAX1 ) 

227  NP=1 

228  DO  1  IN=1 , IRS 

229  CALL  POWERO ( NUE, W,AI ,B0 ,P0) 

230  IF(IN  .EO.  1)  THEN 

231  INIT=C 

23?  ELSE 

233  INIT=1 

mi  END  IF 

235  IF( NP  .  Eq.  NJP )  THEN 

230  GO  TO  909 

237  ELSE 

230  END  IF 
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239  CALL  POWMAI(MDIM,N0E,NMP,INIT,W,PSI,B,AI,IP,NP,L,PMAI) 

240  CALL  P0WPIE(MDIM,NUE,NMP,W,PSI,B,AI,NP,IP,S1,S2,PA,PW) 

241  CALL  JACOBI(MDIM,NUE,NMP ,NJP , W,PSI ,B , IP ,NP ,S1 ,S2,PW,R,Z,P, AI , JC) 

242  CALL  PP0RGM(NUE,W,I,T,AI,B0,q) 

243  CALL  GRADFG(NUE,PO,PA,P,q,G) 

244  CALL  HESIAN(MDIM,NUE,NMP,NJP,W,I,T,PSI,B0,B,IP,NP,P0,PA,S1,S2 

245  A,R,Z,P,q,PW,JC,G,AI,H) 

246  C'M'********************+*************+*****************************'M'**C 

247  C  INITIALIZE  THE  RESIDUAL  C 

248  C****************** ************** ***************** ***************** **  +  *C 

249  DO  2  1=1, NDE 

250  R(I)=-G(I) 

251  2  CONTINUE 

252  IF( IN  .Eq.  1)  THEN 

253  CALL  N0RM22 ( NOE , R , RO ) 

254  ELSE 

255  END  IF 

258  DO  3  1=1, NUE 

257  Z(I)=0.0 

258  3  CONTINUE 

259  C*********** ***********************************+******************#*** *C 

280  C  INITIALIZE  SEARCH  VECTOR  C 

281  C********************************************************+*************C 

282  CALL  MATVEC(MDIM,NUE,H,R,q) 

283  CALL  N0RM22(NUE , Q , QN ) 

284  BE0=1.0/qN 

285  DO  4  1=1,  SUE 

288  P(I)=BEO*q(I) 

287  4  CONTINUE  , 

288  C* ********************************************************* ************C 

289  C  PERFORM  2  CONJUGATE  GRADIENT  ITERATIONS.  C 

270  C***4'******************************************************************C 

271  DO  6  K=1 ,2 

272  C* ************************************************************ *********C 

273  C  UPDATE  PHASE  VECTOR  AND  RESIDUAL  C 

274  C************************ ********************************************* *C 

275  CALL  MATVEC(MDIM,NUE,H,P,q) 

278  CALL  N0RM22(NUE,q,qN) 

277  AK=1 . O/QN 

278  DO  8  1=1, NUE 

279  Z(I)=Z(I)+AK*P(I) 

28..  0  CONTINUE 

20!  DO  7  1=1 , NUE 

2«2  R(I)=R(I)-AK*q(I) 

2«:t  7  CONTINUE 

284  C7*********************************************************************C 

285  C  UPDATE  SEARCH  VECTOR  C 

288  C****************************!*-*****************************************C 
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287  CALL  MATVEC(MDIM,NUE,H,R,q) 

288  CALL  N0RM22 (NUE , Q , QN ) 

289  BEK=1 . O/QH 

290  DO  8  1=1,  NUE 

291  P(I)=P(I)+BEK*q(I) 

292  8  CONTINUE 

293  B  CONTINUE 

294  C*** ************************************************* ******************C 

295  C  KEEP  THE  PHASE  BETWEEN  -PI  AND  PI.  C 

298  C**********************************************************************C 

297  PI= . 31415928E+01 

298  TP= . 628318B3E+01 

299  DO  9  1=1, NUE 

300  W(I)=W(I)+Z(I) 

301  888  CONTINUE 

302  IF(W(I)  .GT.  PI)  THEN 

303  W(I)=W(I)-TP 

304  GO  TO  888 

305  ELSE 

306  END  IF 

307  IF(W(I)  .LT.  -PI)  THEN 

308  W(I)=W(I)+TP 

309  GO  TO  888 

310  ELSE 

311  END  IF 

312  9  CONTINUE 

313  C*********** ************************************************ ***********C 

314  C  COMPUTE  THE  SIDELOBE  LEVEL  AND  COMPARE  WITH  THE  MINIMIUM  ACHIEVED  C 

315  C  IN  ALL  PRIOR  ITERATIONS.  STORE  THE  MINIMUM  IN  A  TEMPORARY  VECTOR.  C 
318  C** ************************************************************ ********C 

317  PDB=10.0*ALOG10(PMAX/P0) 

318  IF(PDB  .LT.  SLLMIN)  THEN 

319  ISLL=IN 

320  SLLMIN=PDB 

321  DO  77  1=1, NUE 

322  TV(I)=W(I) 

323  77  CONTINUE 

324  ELSE 

325  END  IF 

328  CALL  BEAMWD(NUE ,X , Y , AI , W ,HT , HPBW ,BWFN ) 

327  PTL=100.0*(P0/P01) 

32R  1  CONTINUE 

329  POP  CONTINUE 

330  C+*******+**+*+****************************************+****+********+*C 

331  C  CHANGE  THE  COMPUTED  PHASE  FROM  RADIANS  TO  DEGREES  C 

332  C**********************************************************************C 

333  CALL  DEGREE ( NUE , TV , P ) 

33i  Z(1 )=HPBW 
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335  Z(2)=BVFN 

338  Z(3)=ATL 

337  Z(4)=PTL 

338  RETURN 

339  END 

340  C 

341  SUBROUTINE  AMPTAP(NUE, AI ,ATL) 

342  C*+++**+**********+*+***+****+*+**++**++++++********++***++**+**+******C 

343  C  THIS  SUBROUTINE  COMPUTES  THE  AMPLITUDE  TAPER  LOSS.  C 

344  C**  +  *********+4>**********4i*4>*4>*********4>***********  +  ****  +  **  +  *4*********C 

345  REALM  AI(*) 

346  SUM1=0 .0 

347  SUM2=0.0 

348  DO  1  1=1, NUE 

349  SUM1=SUM1+AI(I) 

350  SUM2=SUM2+AI(I)*AI(I) 

361  1  CONTINUE 

352  ATL=100.0*SUM1*SUM1/(NUE*SUM2) 

353  RETURN 

354  END 

355  C 

36«  SUBROUTINE  ANGLES(TMIN,TMAI,PMIN,PMAI,NT,IPM,IPD,U,V) 

357  C***+*****************************+**+****+******************+*********C 
368  C  NT  (14)  -  THE  NUMBER  OF  SAMPLES  IN  THE  THETA  DIRECTION.  C 

359  C  IPM  (14)  -  THE  MINIMUM  NUMBER  OF  SAMPLES  IN  THE  PHI  DIRECTION.  C 

360  C  IPD  (14)  -  THE  INCREMENT  IN  THE  PHI  SAMPLES  FOR  EACH  THETA  VALUE.  C 

361  C****** *4^*********** ft************************************************ +C 

362  C  THE  TOTAL  NUMBER  OF  SAMPLE  POINTS  IS  GIVEN  BT  THE  ARITHMETIC  C 

363  C  PROGRESSION  FORMULA  C 

364  C  C 

385  C  SUM  FROM  K=0  TO  K=NT  OF  (IPM+K*IPD)  EQUALS  NT*(2*IPM+(NT-l)*IPD)/2  C 

386  C****+*****************************************************************C 

367  REALM  U(*)  ,V(*) 

388  RAD= . 17463293E-01 

389  TP= . 628318SE+01 

370  DT=(TMAX-TMIN)/(NT-1) 

371  DR=(SIN(RAD*TMAX)-SIN(RAD*TMIN))/NT 

372  K=0 

373  KI=-IPD 

374  DO  1  1=0, NT-1 

375  KI=KI+IPD 

376  T=TMIN+I*DT 

377  RT=RAD*T 

376  DP=(PMAX-PMIN)/ (KI+IPM-1) 

370  DO  2  J=0, KI+IPM-1 

380  K=K+1 

381  P=PMIN+J*DP 

382  RP=RAD*P 
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383  ST=SIN(RT) 

384  CP=COS(RP) 

385  SP=SIN(RP) 

388  TPS=TP*ST 

387  U(K)=TPS*CP 

388  V(K)=TPS*SP 

389  2  CONTINUE 

390  1  CONTINUE 

391  RETURN 

392  END 

393  C 

394  SUBROUTINE  INITAL(NUE,D ,X ,T , W) 

395  C**********************************************************************C 

398  C  THIS  SUBROUTINE  INITIALIZES  THE  UNKNOWN.  C 

397  C* ************************************************** *******************C 

398  REAL*4  I(*),T(*),W(*) 

399  TP= . 82831863E+01 

400  RAD= . 174B3293E-01 

401  R0=0.6*D 

402  DO  1  1=1 ,NUE 

403  IS=(-1)**I 

404  R=SqRT(I(I)*I(I)+T(I)*T(I)) 

405  W(I)=3.0*RAD*IS 

400  1  CONTINUE 

407  RETURN 

408  END 

409  C 

410  SUBROUTINE  BETASO(HT ,B0) 

4U  o*********************************************************************C 

412  C  THIS  SUBROUTINE  COMPUTES  BETA  AT  THE  MAIN  BEAM  POINT  (UO.VO)  C 

413  C******************* ************************ ***************************C 

414  TP2= . 3B478418E+02 

415  ARGl=HT*SqRT(TP2) 

418  S=SIN(ARG1) 

417  B0=S*S/TP2 

418  RETURN 

419  END 

420  C 

421  SUBROUTINE  BETASI(NMP,U,V,HT,B) 

422  C******************************** ************************************* *C 

423  C  THIS  SUBROUTINE  COMPUTES  BETA  AT  THE  ITH  SAMPLE  POINT  (U(I),V(I))  C 

424  C**** ********* ******************* **************************************0 

■125  REAL*4  U(*),V(*),B(*) 

120  TP2= . 39478418E+02 

127  DO  1  1=1, NMP 

12«  U2=U( I) *U( I) 

129  V2=V(I)*V( I) 

43..  ARGl=HT*SqRT(TP2-U2-V2) 
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431  ARG2=0. 264-0(1) 

433  S=SIN(ARG1 ) 

433  C=C0S(ARG2) 

434  B(I)=S*S*C*C/ (TP2-U2) 

435  1  CONTINUE 

438  RETURN 

437  END 

438  C 

439  SUBROUTINE  PSIMNlCMDIM.NUE.NMP.I.Y.U.V.PSI) 

441  C  THIS  SUBROUTINE  COMPUTES  PSI_{PQ>.  C 

442  C+*****+++**********+****++*+*+++++***++*+*+++*+++*+*+++*+++**+****++++C 

443  REAL*4  PSI(MDIM,NMP),I(*).Y(*),0(*),V(*) 

444  DO  1  1=1, NUE 

445  DO  2  J=1 ,NMP 

44«  PSI(I,J)=COS(X(I)*U( J))*COS(T(I)*V(J)) 

447  2  CONTINUE 

448  1  CONTINUE 

449  RETURN 

450  END 

451  C 

452  SUBROUTINE  POWERO(NUE,W,AI,BO,PO) 

453  C******4^*****************************+********************************C 

454  C  THIS  SUBROUTINE  COMPUTES  THE  MAIN  BEAM  POWER  AT  THE  POINT  (0,0).  C 

458  REAL*4  W(*) ,AI(+) 

457  SUM1=0 .0 

458  SUM2=0 . 0 

469  DO  1  1=1, NUE 

480  SUM1=SUM1+AI(I)*C0S(W(I)) 

481  SUM2=SUM2+AI(I)*SIN(W(I)) 

482  1  CONTINUE 

483  P0=B0*(SUM1*SUM1+SUM2*SUM2) 

484  RETURN 

485  END 
488  C 

487  SUBROUTINE  POWMAI ( MD IM , NUE , NMP , INIT , W , PSI , B , AI , IP , NP , LMAI , PMAI ) 

488  C**+* *****+++* ***************************************+*****************0 

489  C  THIS  SUBROUTINE  COMPUTES  THE  MAXIMUM  POWER  AND  UPDATES  THE  POINTER  C 

470  C  ARRAY  IP  WITH  ALL  THE  PREVIOUS  MAXIMUM  LOCATIONS.  C 

471  C+**  +  ***’l'***’l'**************+****  +  ****'('******  +  ****  +  ****>|t******  +  **>tc*iM<**!tiC 

472  REAL*4  W(*),PSI(MDIM,NMP),B(*),AI(*) 

473  INTEGER*4  IP(*) 

171  C********** ************************************************************£ 
175  C  SORT  THROUGH  THE  DATA  AND  FIND  THE  MAXIMUM  POWER  AND  THE  INDEX  C 

170  C  NUMBER  LMAX.  C 

477  C***********+*******************+**************************************C 

478  PMAX=-1000.0 
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479 

LMAI=-i 

480 

DO  1  1=1, NMP 

481 

SUM1=0.0 

482 

SUM2=0.0 

483 

DO  2  K=1 ,NUE 

464 

SUM1=SUM1+AI(K)*PSI(K,I)*C0S(W(K)) 

485 

SUM2=SUM2+AI(K)+PSI(K,I)+SIN(W(K) ) 

486  2 

CONTINUE 

487 

P=B(I)*(SUM1*SUM1+SUM2*SUM2) 

488 

IF(P  .GT.  PMAX)  THEN 

489 

LMAX=I 

400 

PMAI=P 

491 

ELSE 

492 

END  IF 

493  1 

CONTINUE 

494  C************** *****4>*+**4>**#+*****#*4>****+**‘l'++******  +  ****  +  ********4>**C 

495  C  IF  THIS  IS  THE  FIRST  CALL  THEN  SET  THE  POINTER  TO  THE  FIRST  MAXIMUM  C 

496  C  AND  RETURN.  C 

497  C******** *************** ********************* **************************C 

498  IFflNIT  . EQ.  0)  THEN 

499  IP(1)=LMAI 

500  ELSE 

501  C*********************************************************************+C 


602  C  COMPARE  LMAX  TO  ALL  THE  PREVIOUS  POINTERS  IP(I) ,1=1 . . .NP.  IF  LMAI  C 

503  C  MATCHES  ONE  OF  THE  PREVIOUS  POINTERS  THEN  THE  POINTER  ARRAY  NEED  C 

504  C  NOT  BE  UPDATED.  IF  LMAX  DOES  NOT  MATCH  A  PREVIOUS  POINTER  VALUE  C 

60S  C  THEN  IT  BECOMES  THE  NEIT  VALUE  OF  THE  NP+1TH  ELEMENT  OF  THE  POINTER  C 

soe  C  ARRAY.  C 


507  C**********************************************************************C 

608  IS=0 

609  DO  3  1=1, NP 

510  IF(IP(I)  -Eq.  LMAX)  THEN 

511  IS=1 

612  GO  TO  4 

513  ELSE 

514  END  IF 

sis  3  CONTINUE 

616  4  CONTINUE 

517  IF(IS  .Eq.  0)  THEN 

518  NP=NP+1 

M9  IP(NP)=LMAX 

520  ELSE 

521  END  IF 

522  END  IF 

523  RETURN 

524  END 

525  C 

528  SUBROUTINE  POWPIE( MDIM , NUE , NMP ,W,PSI,B,AI,NP,IP,S1,S2,PA,PW) 
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527  C**************************** **+********************* ******************C 

528  C  THIS  SUBROUTINE  COMPUTES  THE  POWER  AT  THE  ITH  POINT.  C 

529  C* ******* >»******iM<iM<i|i*******************#*****#*****  +  i|i*>Mi**>Mi**********C 

530  REAL*4  W(*) ,PSI(MDIM ,NMP) ,B(*),PW(*),AI(*) 

531  INTEGER*4  IP(*) 

532  S1=0 .0 

533  S2=0.0 

634  DO  1  1  =  1,  NP 

635  INX=IP(I) 

638  SUM1=0.0 

537  SUM2=0.0 

538  DO  2  K=1,NUE 

639  SUMl=SUMi+AI(K)*PSI(K,INI)*COS(W(K)) 

540  SUM2=SUM2+AI(K)*PSI(K, INI)*SIN(W(K)) 

641  2  CONTINUE 

542  PW(I)=B(INX)*(SUM1*SUM1+SUM2*SUM2) 

543  S1=S1+PW(I) 

544  S2=S2+PW(I)*PW(I) 

545  1  CONTINUE 

548  PA=S2/S1 

347  S2=0.G*S2 

548  RETURN 

549  END 

550  C 

551  SUBROUTINE  JACOBI (MDIM , NUE ,NMP ,N JP , W , PSI ,B , IP , NP , SI , S2 ,PW , A .BE 

552  *,PAG,AI,JC) 

553  C**********  +  ***************************************************  +  +  ****,t,*C 

654  C  THIS  SUBROUTINE  COMPUTES  THE  FIRST  DERIVATIVE  OF  PI  WITH  RESPECT  C 

556  C  TO  GAMMA_{M>.  C 

658  C***if +  +  +  +  +  +  +  +  +  +  +  +  + 

557  REAL*4  W(*),PSI(MDIM,NMP),B(*),JC(NJP,MDIM),PW(+),A(*),BE(*) 

668  *,PAG(*),AI(*) 

569  INTEGER*4  IP(«0 

680  Cl=2 .O/Sl 

681  C2=-2.0*S2/(S1*S1) 

582  C******************************** *********** ***************************C 

583  C  COMPUTE  THE  (NP  X  NUE)  JACOBIAN  C 

584  C**********************************************************************C 

586  DO  1  1=1, NP 

588  INX=IP(I) 

587  TB=2 . 0*B( INX) 

588  DO  2  M=1,NUE 

589  SUM1=0 . 0 

570  DO  3  K  =  1 ,  NUE 

571  SUM1=SUM1+AI(K)*PSI(K,INX)*SIN(W(K)-W(M) ) 

572  3  CONTINUE 

573  JC( I , M)=TB*AI(M) *PSI(M , INX) *SUM1 

574  2  CONTINUE 
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575  1  CONTINUE 

57a  c**********************************************************************c 

577  C  SUM  UP  THE  COLUMNS  OF  THE  JACOBIAN  TO  USE  FOR  WORK  VECTORS.  C 

578  C+**++*+*****+***+*+****+**+*++++***********++*+*****++*+*******+****++C 


679 

DO  4  M=1,NUE 

580 

SUM1=0.0 

581 

SUM2=0.0 

582 

DO  B  1=1, NP 

583 

SUM1=SUM1+JC(I,M) 

584 

SUM2=SUM2+PW(I)*JC(I  ,M 

585 

S 

CONTINUE 

580 

A(M)=SUM1 

587 

BE(M)=SUM2 

588 

PAG(M)=C1*SUM2+C2*SUM1 

589 

4 

CONTINUE 

590 

RETURN 

691 

END 

592  C 

593 

SUBROUTINE  GRADFG ( NUE , P 0 

594  C*********** ******** ******************* ********************* ********** *C 

595  C  THIS  SUBROUTINE  COMPUTES  THE  GRADIENT  OF  THE  WEIGHTED  AVERAGE  C 

sea  C  MAIIMUM  SIDELOBE  LEVEL  POWER.  C 

607  C**********************************************************************C 

608  REAL*4  PAG(*) ,POG(+)  ,G(  +  ) 

509  Cl=1.0/P0 

eoo  C2=-PA/(P0*P0) 

eoi  DO  1  M=1 ,NUE 

602  G(M)=C1*PAG(M)+C2*P0G(M) 

803  1  CONTINUE 

604  RETURN 

605  END 


807  SUBROUTINE  PPORGM(NUE,W,X,T,AI,BO,POG) 

808  C^*********************************************************************C 

809  C  THIS  SUBROUTINE  COMPUTES  THE  VECTOR  OF  FIRST  DERIVATIVES  OF  PO  WITH  C 

610  C  RESPECT  TO  GAMMA_{M>  C 

au  c* ************************************************************** *******C 

612  REAL*4  W(*),X(*),T(*),AI(*),POG(*) 

ei3  C=2 . 0*B0 

614  DO  1  M=1 ,NUE 

815  SUM=0 . 0 

8t8  DO  2  K=1 ,NUE 

air  SUM=SUM+AI(K)*SIN(W(K)-W(M) ) 

8i8  2  CONTINUE 

810  POG(M)=C*AI(M)*SUM 

820  1  CONTINUE 

021  RETURN 

822  END 
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fi 


023  C 

624  SUBROUTINE  HESIANCMDIM.NUE.NMP.NJP.W.I.T.PSI.BO.B.IP.NP.PO.PA.Sl 

625  *,S2,A,BE,PAG,P0G,PW,JC,G,AI>H) 

626  C**********************************************************************^ 

627  C  THIS  SUBROUTINE  COMPUTES  THE  HESSIAN  OF  THE  WEIGHTED  AVERAGE  C 

628  C  MAXIMUM  SIDELOBE  LEVEL  POWER.  C 

628  C**********************************************************************C 

630  REAL*4  H(MDIM,MDIM).W(*),X(*),T(*),PSI(MDIM,NMP),B(*) 

631  AIA(*),BE(*),PAG(*),POG(*),PW(*),JC(NJP,MDIM).G(*),AI(*) 

632  INTEGER+4  IP(*) 

633  C1=-2.0*P0/(S1*S1) 

634  C2=-l . 0/ (PO*PO) 

635  C3=2.0*B0 

837  C  COMPUTE  ALL  INFORMATION  FOR  THE  MTH  ROW.  C 

838  C't'**  +  ****4"**>f***4'>|i**'t>**4>***it>***4i4>*********4>4i***#*>f***4i*H<  +  *****4'*#****  +  #C 

638  DO  2  M=1 ,NUE 

640  CM=A(M)+PAG(M) 

842  C  COMPUTE  ALL  INFORMATION  FOR  THE  NTH  COLUMN.  C 

643  C***********4^************ ******************************************** *C 


644  DO  3  N=M,NUE 

845  DWNM=W(N)-W(M) 

848  AIMN=AI(M)+AI(N) 

847  CDWNM=COS (DWNM) 

848  P0MN=C3*AIMN*CDWNM 

848  P2=PA*P0MN 

860  SUM1=0 .0 

861  SUM2=0.0 

652  DO  4  1=1, NP 

853  INX=IP(I) 

654  PIMN=2 . 0*B  ( INX)  ♦PSI(N ,  INX)  *PSI  (M ,  INX)*CDWNM 

655  SUM1=SUM1+JC(I,M)*JC(I,N)+PW(I)*PIMN 

656  SUM2=SUM2+PIMN 

857  4  CONTINUE 

858  P1=P0*(G(N)*P0G(M)+G(M)*P0G(N)) 

658  P3=C1+(A(N)*BE(M)-CM*BE(N)+S1*SUM1-S2+AIMN*SUM2) 

660  H(M,N)=C2*(P1+P2+P3) 

661  3  CONTINUE 

882  2  CONTINUE 

003  C*************************************************+***********+********C 

604  C  COMPLETE  THE  BOTTOM  TRIANGLE  OF  THE  SYMMETRIC  MATRII .  C 

065  C**********************************************************************C 
060  DO  6  M=2 ,NUE 

007  DO  7  N=1  ,M-1 

006  H(M , N) =H(N ,M) 

069  7  CONTINUE 

07o  fl  CONTINUE 
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fl7i  RETURN 

072  END 

673  C 

874  SUBROUTINE  N0RM22(N,A, AN) 

675  C********************  ***************************************  *****4"H"t|**C 

678  C  THIS  SUBROUTINE  COMPUTES  THE  EUCLIDIAN  NORM  SQUARED  OF  A.  C 

677  C**********************************************************************C 
878  REAL*4  A( * ) 

679  AN=0.0 

680  DO  1  1=1 ,N 

681  AN=AN+A(I)*A(I) 

882  1  CONTINUE 

683  RETURN 

684  END 

885  C 

886  SUBROUTINE  MATVEC(MDIM,NUE,H,P,Q) 

687  C*************************************** ***************************** **C 
888  C  THIS  SUBROUTINE  COMPUTES  THE  EUCLIDIAN  NORM  SQUARED  OF  A.  C 

689  C**** ***************************************************** *************C 
890  REAL*4  H(MDIM,MDIM),P(*),Q(*) 

691  DO  1  M=1 .HUE 

692  SUM=0 .0 

693  DO  2  N=1,NUE 

894  SUM=SUM+H(M ,N)*P(N) 

895  2  CONTINUE 

898  q(M)=SUM 

897  1  CONTINUE 

898  RETURN 

899  END 

700  C 

701  SUBROUTINE  DEGREE (NUE,W,R) 

702  C**********************************************************************C 

703  C  THIS  SUBROUTINE  CONVERTS  RADIANS  INTO  DEGREES  C 

704  C**********************************************************************C 

706  REAL*4  W(*),R(*) 

708  PI= . 3141SD20E+O1 

707  DO  1  1=1, NUE 

705  R(I)=(180.0/PI)*W(I) 

709  1  CONTINUE 

710  DO  2  1=1, NUE 

711  3  CONTINUE 

712  IF(R(I)  .LT.  -180.0)  THEN 

713  R( I )=R( I )+380 . 0 

714  GO  TO  3 

715  ELSE  IF(R(I)  .GT.  180.0)  THEN 
R( I)=R( I ) -380 . 0 
GO  TO  3 

ELSE 


718 

717 

718 
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718 

END  IF 

720  2 

CONTINUE 

721 

RETURN 

722 

END 

723  C 

724 

SUBROUTINE  BEAMWD ( NUE , I , T , AI , W , HT , HPBW , BWFN) 

725  C*’t‘  +  ’t,*************  +  *****,*‘**t,*  +  ,tl**  +  *****  +  ****  +  ***l,,*,,,************','***',",",",,C 

726  C  THIS  SUBROUTINE  COMPUTES  THE  HALF  POWER  BEAMWIDTH  AND  THE  BBAMWIDTH  C 

727  C  BETWEEN  FIRST  NULLS  OF  THE  POWER  PATTERN.  C 

7j8  C***************************’*******************************************C 

728  REAL*4  I(*),T(*),AI(*),W(*) 

730  RAD= . 174E3293E-01 

731  TP= . 628318B3E+01 

732  TP2= . 39478418E+02 

733  C*+*++***+********+++*+****+**+*****+**+***+*+*********+***************C 

734  C  COMPUTE  THE  POWER  IN  THE  MAIN  BEAM.  C 

735  C**********************************************************************C 

738  ARGi=HT*TP 

737  S=SIN(ARG1) 

738  B0=S*S/TP2 

739  SUM1=0.0 

740  SUM2=0 .0 

741  DO  4  1=1 ,NUE 

742  SUM1=SUM1+AI(I)+C0S(W(I)) 

743  SUM2=SUM2+AI(I)*SIN(W(I)) 

744  4  CONTINUE 

745  P0=B0+(SUM1*SUM1+SUM2*SUM2) 

748  C*******4'****4'********4|*  +  +  **********4'II'*************I',*’(,**I,,**I,,********',,**C 

747  C  FIND  THE  HALF  POWER  POINT  AND  THE  FIRST  NULL  POWER  POINT.  C 

748  C*+****+**+++++++*******+++++****++++***+****+*****+*+++**+*+**********C 


740 

T=0 . 0 

750 

DT= . OB 

751 

PL0WER=1 . 0 

752 

ITAG=0 

753 

P=46 . 0 

754 

RP=RAD*P 

756 

TPCP=TP*COS (RP ) 

750 

TPSP=TP*SIN(RP) 

757  5 

CONTINUE 

768 

T=T+DT 

759 

R=TP*SIN(RAD*T ) 

760 

ST=SIN(RAD*T) 

701 

U=ST+TPCP 

702 

V=ST*TPSP 

703 

U2=U*U 

704 

V2=V*V 

705 

ARG1=HT*SQRT(TP2-U2-V2) 

700 

ARG2=0 . 26*U 
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767  S=SIN(ARG1) 

768  C=C0S(ARG2) 

769  BI=S+S*OC/ (TP2-U2) 

770  S1=0.0 

771  S2=0 . 0 

772  DO  6  K=1 ,NUE 

773  UI=U*X(K) 

774  7T=V+T(K) 

775  PSI=COS(X(K)*U)*COS(T(K)*V) 

776  S1=S1+AI(K)*PSI*C0S(W(K) ) 

777  S2=S2+AI(K)*PSI*SIN(W(K)) 

778  6  CONTINUE 

779  P0I=BI*(S1*S1+S2*S2)/P0 

780  IF(POI  .LE.  PLOWER)  THEN 

781  IF( (POI  .LE.  0 . B ) . AND . ( ITAG  .EQ.  0))  THEN 

782  HPBW=2.0*T 

783  ITAG=1 

784  ELSE 

785  END  IF 

788  PLOWER=POI 

787  GO  TO  B 

788  ELSE 

789  BWFN=2 . 0*T 

790  GO  TO  7 

791  END  IF 

792  7  CONTINUE 

793  RETURN 

794  END 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-field-of-view  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


